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Abstract 

It is shown that max-stable random vectors in [0, 00)°' with unit Frechet marginals 
are in one to one correspondence with convex sets K in [0, oo)'^ called max-zonoids. 
The max-zonoids can be characterised as sets obtained as limits of Minkowski sums 
of cross-polytopes or, alternatively, as the selection expectation of a random cross- 
polytope whose distribution is controlled by the spectral measure of the max-stable 
random vector. Furthermore, the cumulative distribution function P{^ < x} of a 
max-stable random vector ^ with unit Frechet marginals is determined by the norm 
of the inverse to x, where all possible norms are given by the support functions of 
(normalised) max-zonoids. As an application, geometrical interpretations of a number 
of well-known concepts from the theory of multivariate extreme values and copulas are 
provided. 

Keywords: copula; max-stable random vector; norm; cross-polytope; spectral measure; 
support function; zonoid 
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1 Introduction 



A random vector ^ in M'^ is said to have a max-stable distribution if, for every n > 2, the 
cooordinatewise maximum of n i.i.d. copies of ^ coincides in distribution with an affine 
transform of ^, i.e. 

^(l)V--Ve(")~ane + &n (1-1) 

for a„ > and bn G M'^. If (11. ip holds with 6„ = for all n, then ^ is called strictly 
max-stable, see, e.g., P [22| [28]. 

Since every max-stable random vector ^ is infinitely divisible with respect to coordinate- 
wise maximum, its cumulative distribution function satisfies 

= P{^ < x} = <^ ' J^i' - ' xGM^ (1.2) 

I , otherwise , 

where a G [— oo, oo)'^, the superscript c denotes the complement and /i is a measure on 
[a, oo] \ {a} called the exponent measure of ^, see [281 Prop. 5.8]. Note that all inequalities 
and segments (intervals) for vectors are understood coordinatewise. 

Representation (11.21) shows that the cumulative distribution function of ^ can be repre- 
sented as the exponential F(x) = e~'^^^^ of another function u. If ^ is strictly max-stable 
and a = 0, then u is homogeneous, i.e. z/(sx) = s~"i'{x) for all s > and some a > 0. This 
fact can be also derived from general results concerning semigroup- valued random elements 
[6]. If a = 1, an example of such function i>{x) is provided by z/(x) = ||x*||, i.e. a norm of 
X* = (xi^, . . . , x^^) for X = (xi, . . . , Xd) G [0, oo)'^. One of the main aims of this paper is to 
show that this is the only possibility and to characterise all norms that give rise to strictly 
max-stable distributions with a = 1. 

Every norm is homogeneous and sublinear. It is known [331 Th. 1.7.1] that each bounded 
homogeneous and sublinear function g on M."^ can be described as the support function of a 
certain convex compact set K, i.e. 

g{x) = h{K,x) = snp{{x,y) : y E K} , 

where {x,y) is the scalar product of x and y. In Section [2] we show that every standardised 
strictly max-stable distribution with a = 1 is associated with the unique compact convex 
set K C [0, oo)'^ called the dependency set. The dependency sets are suitably rescaled sets 
from the family of sets called max-zonoids. While classical zonoids appear as limits for 
the sums of segments [SHI Sec. 3.5], max-zonoids are limits of the sums of cross-polytopes. 
The contributions of particular cross-polytopes to this sum are controlled by the spectral 
measure of the max-stable random vector. It is shown that not every convex compact set 
for d > 3 corresponds to a strictly max-stable distribution, while ii d = 2, then the family 
of dependency sets is the family of all "standardised" convex sets, see also [TT] for the 
treatment of the bivariate case. This, in particular, shows a substantial difference between 
possible dependency structures for bivariate extremes on one hand and multivariate extremes 
in dimensions three and more on the other hand. 
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The geometrical interpretation of max-stable distributions opens a possibility to use tools 
from convex geometry in the framework of the theory of extreme values. For instance, the 
polar sets to the dependency set K appear as multivariate quantiles of the corresponding 
max-stable random vector, i.e. the level sets of its cumulative distribution function. In 
the other direction, some useful families of extreme values distributions may be used to 
construct new norms in which acquire an explicit probabilistic interpretation. The norms 
corresponding to max-stable distributions are considered in Section [31 

Section H] deals with relationships between spectral measures of max-stable laws and 
geometric properties of the corresponding dependency set. In Section [5] it is shown that 
a number of dependency concepts for max-stable random vectors can be expressed using 
geometric functionals of the dependency set and its polar. Here also relationships to copulas 
are considered. It is shown that max-zonoids are only those convex sets whose support 
functions generate multivariate extreme value copulas. 

It is well known that Z is (classical) zonoid if and only if e~'*(^'^'^ is positive definite, 
see [331 P- 194]. In Section [6] we establish a similar result for the positive definiteness of 
the exponential with respect to the coordinatewise maximum operation in case Z is a max- 
zonoid. 

Section[7]describes some relationships between operations with convex sets and operations 
with max-stable random vectors. Finally, Section [8] briefly mentions an infinite-dimensional 
extension for max-stable sample continuous random functions. 



2 Dependency sets and max-zonoids 

Let ^ be a max-stable random vector with non-degenerate marginals. By an affine transfor- 
mation it is possible to standardise the marginals of ^, so that ^ has $q (Frechet distributed) 
marginals, where 

^ , , 1 0, X < , 

^ ^ |e-^ , x>0, 

or (WeibuU or negative exponential distributed) marginals, i.e. 

^„ X = < a>0, 
\l, a;>0, 

or A (Gumbel or double exponentially distributed) marginals, i.e. 

A(x) = exp{-e"''} , x G M . 

By using (possibly non-linear) monotonic transformations applied to the individual coordi- 
nates it is possible to assume that all marginals are $i, see ^281 Prop. 5.10] and [21 Sec. 8.2.2]. 
In this case we say that ^ has unit Frechet marginals or has a simple max-stable distribu- 
tion, see also [lO]. Sometimes we say that ^ = (^i,...,^d) has a semi-simple max-stable 
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distribution if its rescaled version (ci^i, . . . yCdCd) has a simple max-stable distribution for 
some Ci, . . . , Q > 0. 

If ^ has a simple max-stable distribution, then P2i Prop. 5.11] implies that the exponent 
in (11.21) has the following representation 

z/(x) =/i([0,x]^) = / max I —] a{da) , x G [0, cx)]'^ \ {0} , (2.1) 

Jg^ l<i<d \XiJ 

where 

[0,x] = xf^jO,Xi] , X = (xi, . . . ,Xd) , 

S+ = {s G E : llsll = 1} is a sphere in E = [0, oo)'^ with respect to any chosen norm (from 
now on called the reference sphere and the reference norm) and a is a finite measure on §+ 
(called the spectral measure of ^) such that 

/ aio{da) = l, i = l,...,d. (2.2) 

A similar representation is described in Th. 4.3.1] for the special case of S+ being the 
unit simplex. 

We now aim to relate the function v{x) from (12.11) to the support function of a certain 
compact convex set. Recall that the support function of a set M C is defined as 

h{M, x) = sup{(2, x) : ze M}, 

where {z, x) is the scalar product in W^. Let ei, . . . , be the standard orthonormal basis in 
W^. For every a = (ai, . . . , Od) G consider the cross-polytope 

Aa = conv({0, aid, • • • , adCd}) , 

where conv(-) denotes the convex hull of the corresponding set. Note that conv({aiei, . . . , adCd}) 
is a simplex. Then 

h{Aa,x) = h{A^,a) = max (a^Xj) 

l<i<d 

for every a G S+ and x G E. For x = (xi, . . . , Xd) G E write x* = (xj~^, . . . , x^^). Then (12.11) 
can be expressed as 

iy{x*)= / h{Aa, x)a{da) , xgE. (2.3) 

The function /(x) = z/(x*) is called the stable tail dependence function, see [21 p. 257]. 

It is well known that the arithmetic sum of support functions of two convex compact sets 
K and L is the support function of their Minkowski sum 

K + L = {x + y: x e K, y e L} , 

i.e. h{K+L, x) equals h{K, x) + h{L, x). Extending this idea to integrals of support functions 
leads to the expectation concept for random convex compact sets, see [1] and [231 Sec. 2.1]. If 
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X is a random compact convex set [23] such that ||X|| = sup{||x|| : x G X} is integrable, then 
the selection expectation (also called the Aumann expectation) of X is the set of expectations 
of for all random vectors ^ such that ^ G X a.s. If the underlying probability space is 
non-atomic, or X is a.s. convex, then EX is the unique compact convex set that satisfies 

E/i(X,x) = /i(EX,x) 

for all X, see [231 Th. 11.1.22]. 

Let 0"! be the spectral measure a normalised to have the total mass 1. If is distributed 
on §+ according to ui, then is a random convex compact set whose selection expectation 
satisfies 

/i(E Z\^, x) = [ h{Aa, x)a{da) . (2.4) 

Condition (12.21) further implies that 

a(§+)/i(EZ\^,ei) = 1, i = l,...,d. (2.5) 

Since /i(EZ\^,ej) = E/i(Z\^,ej) = Er^j, we have a(E>+)Erii = 1 for i = l,...,d. Together 
with (12.1 p and (II. 2p these reasons lead to the following result. 

Theorem 2.1. A random vector ^ is max-stable with unit Frechet marginals if and only if 
its cumulative distribution function F{x) = P{^ < x} satisfies 

F{x) = exp{-ch{EAr„x*)} , xgE, 

for a constant c > and a random vector G S+ such that cErj = (1,...,1). 
U K = cEAr„ then 

F(a:) = e-'^^^'"*) , xeE. (2.6) 
Furthermore, note that K = E Acri with E(c?7) = (1, . . . , 1). 

Definition 2.2. The set K = cEA^ where c > and r/ is a random vector on S+ is said to 
be a max-zonoid. If (Xi is the distribution of r], then a = ca\ is the spectral measure of K. If 
cErj = (1, . . . , 1), then the max-zonoid K is called the dependency set associated with the 
spectral measure a (or associated with the corresponding simple max-stable random vector). 

Proposition 2.3. A convex set K is a max-zonoid if and only if there exists a semi-simple 
max-stable vector ^ with cumulative distribution function F{x) = e~^^^'^*^ for all s G E. 

Proof. Sufficiency A semi-simple max-stable ^ can be obtained as ^ = a^' = {ai^[, . . . , adCd) 
for simple max-stable vector ^' and a = (ai, . . . , ad) G (0, oo)'^. Let K' be the dependency 
set of By Theorem 12.11 

P{^ <x} = P{a^' <x} = e-'^^^''""*) = e-^(^'"*) , x G E , 

for K = aK' = {(fliXi, . . . , adXd) : (xi, . . . , x^) G K}. 
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Necessity. If is a max-zonoid, then K' = aK is a dependency set for some a G (0, 00)"^. 
If ^' is max-stable with dependency set K', then it is easily seen that a^' has the cumulative 
distribution function □ 

Proposition 12.31 means that each max-zonoid can be rescaled to become a dependency 

set. 

Proposition 2.4. A max-zonoid K always satisfies 

A,cKc[0,z] (2.7) 

for some 2; G E. 

Proof. The result follows from the following bound on the support function of E 

d 

h{Ay, x) = max 'E{riiXi) < E h{A^, a;) < E N rjiXi = h{[0, y],x) 

l<i<d ^ 

2=1 

where y = Er], so that (12. 7p holds with z = cy. □ 

The normalisation condition (12. 5p and (12. 7p imply that the dependency set of a simple 
max-stable distribution satisfies 

Z\(i,...,i) = conv{0, ei, . . . , e,} C C [0, l]'^ , (2.8) 

where is called the unit cross-polytope. 

The selection expectation of has the support function given by 

h{EArj,x)= \\{aixi, . . . ,adXd)\\oc,<y{da) , (2.9) 

where || ■ ||oo is the £oo-norm in M'^. If the £oo-iiorm in (12.90 is replaced by the £i-norm, i.e. 
the absolute value of the sum of the coordinates and integration is carried over the whole 
sphere, then (12.91) yields the support function of a zonoid, see [321 Sec. 3.5]. This provides 
one of the reasons for calling E a max-zonoid. Note that max-zonoids form a sub-family 
of sets called (i-zonoids in [29] . 

It is possible to define a max-zonoid as the selection expectation of where ( is any 
random vector in E (not necessarily on S_|_). The corresponding spectral measure a on 
can be found from 

g{a)a{da) = E [ilCII^(^)] (2.10) 
for all integrable functions (7 on S+. Indeed, 

h{A^,x)a{du) = E[||C||M^c/IICII>^)] = ^HA^,x) . 



6 



If all coordinates of ( have the unit mean, then the selection expectation of becomes a 
dependency set. 

An alternative representation of max-stable laws [28, Prop. 5.11] yields that 

F{x) = exp |- max (^^, . . . , (2.11) 

for non-negative integrable functions fi, . . . , fa satisfying 

/ fi{s)ds = l, i = l,...,d. 
Jo 

Thus ^ 

h{K,x)= / max{fi{s)xi,...,fd{s)xd)ds, 
Jo 

i.e. the dependency set K is given by the selection expectation of the cross-polytope 
where /(r/) = (/i(?7), . . . , fdiv)) r] is uniformly distributed on [0, 1]. The corresponding 
spectral measure can be found from fl2.10p for ( = firf). 

Theorem 2.5. li d = 2, then each convex set K satisfying Ii2.8\) is the dependency set of 
a simple max-stable distribution. If d > 3, then only those K that satisfy and are 
max-zonoids correspond to simple max-stable distributions. 

Proof. Consider a planar convex polygon K satisfying 02.81) . so that its vertices are a° = 
Ci, a^, . . . , a™ = 62 in the anticlockwise order. Then K equals the sum of triangles with 
vertices (0, 0), {a\~^ — a\, 0), (0, — al~^) for i = 1, . . . , m, where a* = {a\, al). Thus (12.31) 
holds with a having atoms at with mass where Ui = {a\~^ — a\,a2 — ^2"^) 

i = 1, . . . ,m. The approximation by polytopes yields that a general convex K satisfying 
(12. 8p can be represented as the expectation of a random cross-polytope and so corresponds 
to a simple max-stable distribution. 

Theorem 12.11 implies that all max-zonoids satisfying (12.81) correspond to simple max- 
stable distributions. It remains to show that not every convex set K satisfying (12.81) is a 
dependency set in dimension d > 3. For instance, consider set L in which is the convex 
hull of 0,61,62,63 and (2/3,2/3,2/3). All its 2-dimensional faces are triangles, so that this 
set is indecomposable by [T71 Th. 15.3]. Since L is a polytope, but not a cross-polytope, it 
cannot be represented as a sum of cross-polytopes and so is not a max-zonoid. □ 

The support function of the dependency set K equals the tail dependence function (12. 3p . 
If an estimate /(■) of the tail dependence function is given for a finite set of directions 
Ui, . . . ,Um, it is possible to estimate K, e.g. as the intersection of half-spaces {x e E : 
{x,ai) < /(cti)}. However, this estimate should be use very cautiously, since the obtained 
polytope K is not necessarily a max-zonoid in dimensions three and more. While this 
approach is justified in the bivariate case (see also [H]), in general, it is better to use an 
estimate a of the spectral measure a in order to come up with an estimator of K as 

h{K,x)= / h{Aa, x)a{da) . 
Js+ 
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Being the expectation of a cross-polytope, the obtained set is necessarily a max-zonoid. 
The set 

K° = {xEE: h{K,x) < 1} 

is called the polar (or dual) set to K in E, see [331 Sec. 1.6] for the conventional definition 
where E is replaced by M'^. If K is convex and satisfies (12. 8p . then its polar K" is also convex 
and satisfies the same condition. Furthermore, 

{x e E : F{x) > a} = {x G E : e^^^^'^*) > a} 

= {x* : X eE, h{K,x) < - log a} 
= (-loga){a;* : x E K°} , 

i.e. multivariate quantiles of the cumulative distribution function of a simple max-stable 
random vector are inverted rescaled variants of the polar set to the dependency set K. The 
level sets of multivariate extreme values distributions have been studied in |8J. Note that 
the dimension effect described in Theorem 12.51 restricts the family of sets that might appear 
as multivariate quantiles in dimensions d > 3. 

The ordering of dependency sets by inclusion corresponds to the stochastic ordering of 
simple max-stable random vectors, i.e. if ^' and ^" have dependency sets K' and K" with 
K' C K", then P{^' < x} > P{C < x} for all x G E. 

A metric on the family of dependency sets may be used to measure the distance between 
random vectors ^' and ^' with simple max-stable distributions. Such distance can be defined 
as the Hausdorff distance between the dependency sets of ^ and ^' or any other metric 
for convex sets (e.g. the Lebesgue measure of the symmetric difference or the Lp-distance 
between the support functions). In the spirit of the Banach-Mazur metric for convex sets 
(or linear spaces), a distance between two dependency sets K' and K" can be defined as 

d 

m{K', K") = loginf{ JJ A, : K' G XK", K" C XK', X G (0, oo)'^} , 
1=1 

where XK = {(AiXi, . . . , A^x^) : (xi, . . . , x^) G K} with A = (Ai, . . . , A^). If ^' and C," have 
dependency sets K' and K" respectively, then m{K\ K") is the logarithm of the smallest 
value of (Ai ■ ■ ■ A^) such that ^' is stochastically smaller than A^" and ^" is stochastically 
smaller than A^'. For instance, the distance between the unit cross-polytope and the unit 
square (for ci = 2) is log 4, which is the largest possible distance between two simple bivariate 
max-stable laws. 

3 Norms associated with max-stable distributions 

Note that the support function of a compact set L is sublinear, i.e. it is homogeneous and 
subadditive. If L is convex symmetric and contains the origin in its interior, then its support 
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function h{L, x) defines a norm in R'^. Conversely, every norm defines a symmetric convex 
compact set in M.^ with the origin in its interior, see Th. 15.2]. 

Let K he a convex set satisfying fl2.8l) . The corresponding norm || ■ ||x can be defined as 
the support function of the set L obtained as the union of all symmetries of K with respect 
to coordinate planes, i.e. 

\\x\\k = h{L,x) = h{K,\x\) , xeM.'^, 

where |x| = (|a;i|, . . . , The norm \\x\\k is said to be generated by the max-zonoid 

K. Note that the origin belongs to the interior of L and \\x\\k = h{K,x) for x e E. The 
following result shows that distributions of max-stable vectors correspond to norms generated 
by max-zonoids. 

Theorem 3.1. Let \\ ■ \\ be a norm on M."^. The function 

= exp{-||x*||}, xeE, (3.1) 

is the cumulative distribution function of a random vector ^ in E if and only if \\x\\ = h{K, \x\) 
is the norm generated by a max-zonoid K. In this case the random vector ^ is necessarily 
semi-simple max-stable. 

Proof. Sufficiency. If ii" is a max-zonoid. Proposition 12.31 implies that there exists a semi- 
simple max-stable vector ^ with cumulative distribution function e~^^^'^*^ = e~"^*". 
Necessity. If (13. ip is the cumulative distribution function of a random vector ^, then 

P{e^^^ V ■ ■ ■ V e^") <x} = e-"H^*ll = e-ll('^"'^)*ll = < n-^x} 

for all X G E and i.i.d. copies . . . , of ^. Thus, ^ is necessarily semi-simple max-stable. 
Proposition 12.31 implies that (13. ip holds with the norm generated by a max-zonoid K. □ 

The space M.'^ with the norm || ■ \\k becomes a finite-dimensional normed linear space, 
also called the Minkowski space, see [36]. If this is an inner product space, then the norm 
is necessarily Euclidean. Indeed, if K is the intersection of a centred ellipsoid with E and 
satisfies (12.80 . then this ellipsoid is necessarily the unit ball. 

Another common way to standardise the marginals of a multivariate extreme value distri- 
bution is to bring them to the reverse exponential distribution (or unit WeibuU distribution) , 
see [121 Sec. 4.1]. In this case, the cumulative distribution function turns out to be 

F(x) = e-ll^ll^ xG(-oo,0]'^. 

The fact that every max-stable distribution with reverse exponential marginals gives rise to 
a norm has been noticed in [121 P- 127], however without giving a characterisation of these 
norms. 

Note that the norm of \\x\\x can be expressed as 

||x||a' = ||a;|| \\u:^\\k , 
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where = belongs to the reference sphere §+. If the reference norm is ii, then §_|_ 

is the unit simplex and the norm ||-U3;||;^ of m = (ti, . . . , t^^i, 1 — ti — ■ ■ ■ — td-i) G S+ can be 
represented as a function A(ti, . . . , td-i) of ti, . . . , t^-i > such that ti + ■ ■ ■ + t^-i < 1. If 
d = 2, then A{t), < t < 1, is called the Pickands function, see |22] and for the multivariate 
case also [I3ll22]. In general, the norm u G §+, is an analogue of the Pickands function. 

Example 3.2. The dependency set K being the unit cube [0, l]'^ (so that \\x\\k is the £i-norm) 
corresponds to the independence case, i.e. independent coordinates of ^ = (^i, . . . ,^d)- The 
corresponding spectral measure allocates unit atoms to the points from the coordinate axes. 

Furthermore, K being the unit cross-polytope (so that \\x\\k is the £oo-norm) gives rise 
to the random vector ^ = (^i, . . . with all identical $i-distributed coordinates, i.e. the 
completely dependent random vector. The corresponding spectral measure has its only atom 
at the point from §+ having all equal coordinates. Note that the unit cube is dual (or polar) 
set to the unit cross-polytope. 

Example 3.3. The ip-noim \\x\\p with p > 1 generates the symmetric logistic distribution [21 
(9.11)] with parameter a = 1/p. The strength of dependency increases with p. 

Example 3.4. A useful family of simple max-stable bivariate distribution appears if the 
functions /i, /2 in 02.111) are chosen to be the density functions of normal distributions, see 
[22| Sec. 3.4.5] and [2^ p. 309]. It is shown in [20] that these distributions appear as limiting 
distributions for maxima of bivariate i.i.d. Gaussian random vectors. The corresponding 
norm (which we call the Hiisler-Reiss norm) is given by 

\\x\\k = Xi<^{\ + ^ log — ) + X2^\ - ^ log — ) , 

2A X2 2A X2 

where A G [0, cxd]. The cases A = and X = oo correspond to complete dependence and 
independence, respectively. 

Example 3.5. The bivariate symmetric negative logistic distribution ^ p. 307] corresponds 
to the norm given by 

\\x\\k = — M\x\\p , 

where A G [0, 1] and p G [— oo, 0]. 



4 Spectral measures 

Since the dependency set determines uniquely the distribution of a simple max-stable random 
vector, there is one to one correspondence between dependency sets and normalised spectral 
measures. It is possible to extend this correspondence to max-zonoids on one side and all 
finite measures on S+ on the other one, since both uniquely identify semi-simple max-stable 
distributions. While the spectral measure depends on the choice of the reference norm, the 
dependency set remains the same whatever the reference norm is. 

It is shown in [5] that the densities of the spectral measure on the reference simplex 
S+ = convjei, . . . , e^} can be obtained by differentiating the function z/(a;) = fi{[0,xY) for 
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the exponent measure /i. This is possible if the spectral measure is absolutely continuous 
with respect to the surface area measures on relative interiors of all faces of the simplex and, 
possibly, has atoms at the vertices of S+. Following the proof of this fact in [21 Sec. 8.6.1], 
we see that 

lim DAiy{z) = (-l)l^l"^L'A/i({a: G E : xj > zjj G A; xj = 0, j ^ A}) , 

where A C {1, . . . , (i}, |A| is the cardinality of A, and denotes the mixed partial derivative 
with respect to the coordinates with numbers from A. 

The derivatives of u can be expressed by means of the derivatives of the stable tail 
dependence function l{z) = u^z*) = \\z\\k- Indeed, 

DAu{z)=DAl{z*)i-l)^^^l[zj\ 

Thus, the densities of the exponent measure can be found from 
DAfi{{x G E : Xj < Zj,j G A; Xj = 0,j ^ A}) = (-l)!^!"^ lim DAliz*)Y[zj'^ . 

In particular, the density of in the interior of E can be found from the rfth mixed partial 
derivative of the norm as 

After decomposing these densities into the radial and directional parts, it is possible to obtain 
the spectral measure by 

a{G) = fi{{tu : ueG,t>l})= [ f{z)dz 

J{tu: MGG,t>l} 

for every measurable G from the relative interior of S+. The relationship between spectral 
measures on two different reference spheres is given in [2, p. 264]. 

Proposition 4.1. A d-times continuously differentiable function l{x), x & E, is the tail 
dependency function of a simple niax-stable distribution if and only if I is sublinear, takes 
value 1 on all basis vectors, and all its mixed derivatives of even orders are non-positive and 
of odd orders are non-negative. 

Proof. The necessity follows from Theorem 12.11 and the non-negativity condition on the 
exponent measure /i. In the other direction, the sublinearity property implies that / is the 
support function of a certain convex set K, see [33| Th. 1.7.1]. The condition on the sign of 
mixed derivatives yields that the corresponding densities of /i are non-negative, i.e. K is the 
max-zonoid corresponding to a certain spectral measure. □ 
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In the planar case, [33l Th. 1.7.2] implies that the second mixed derivative of the (smooth) 
support function is always non-positive. Accordingly, all smooth planar convex sets satisfying 
fl2.8p are dependency sets. 

A number of interesting measures on the unit sphere appear as curvature measures of 
convex sets [331 Sec. 4.2]. A complete interpretation of these curvature measures is possible 
in the planar case, where the curvature measure becomes the length measure. The length 
measure Si (L, A) generated by a smooth set L associates with every measurable A C §^ 
the 1-dimensional Hausdorff measure of the boundary of L with unit normals from A. The 
length measure for a general L is defined by approximation. Recall that is the part of 
the unit circle lying in the first quadrant. 

Theorem 4.2. A measure a on E>\. is the spectral measure of a simple max-stable random 
vector ^ with dependency set K if and only if a is the restriction on S>]_ of the length measure 
generated by K = {(xi,X2) : {x2,Xi) G K} with K satisfying 

Proof. Sufficiency. Consider a planar convex set K satisfying (12.81) . Let a{da) be the length 
measure of L = K, i.e. the first-order curvature measure Si{L,da). Then 

/ hiAa,x)a{da)= [ h{A,,a)Si{L,da) = 2ViA,,L) , 

where V{Ax,L) denotes the mixed volume (the mixed area in the planar case) of the sets 
A^ and L, i.e. 

2ViA,, L) = V2iL + A^) - V2iL) - V2iA^) , 

see p3| Sec. 5.1]. Because of (12.81) . the integral over the full circle with respect to cr 
coincides with the integral over S-*^ fl [0, oo)^ = Sj*^. It remains to show that 2V{Ax,L) 
equals h{K,x). If z = (-21,-22) G L is any support point of L in direction x = {xi,X2), i.e. 
h{L, x) = {z, x), then 

2V{A^, L) = Z1X2 + Z2X1 = h{K, x) . 

An alternative proof follows the construction from Theorem 12. 5[ Indeed, a polygonal K 
can be obtained as the sum of triangles. A triangle -^(4,5) with vertices (0, 0), (t, 0) and (0, s) 
corresponds to the spectral measure having the atom at (t, s)c~^ with mass c = \^W+~s^ 
and therefore coincides with the length measure of ^(s,t) = ^it,s) restricted onto §J',_. Since 
the spectral measure of K is the sum of spectral measures of these triangles, it can be 
alternatively represented as the sum of the length measures. A general K can be then 
approximated by polygons. 

Necessity. Assume that a measure a on §^ is the spectral measure of a simple max-stable 
law with dependency set K. If now a' is chosen to be the length measure of K restricted 
onto S^, then a' generates the max-zonoid K. Finally, a = a' by the uniqueness of the 
spectral measure. □ 

Theorem 14.21 implies that the length of the boundary of K inside (0, oo)"^ equals the total 
mass of the spectral measure on Given (12.81) . an obvious bound on this boundary length 
implies that this total mass lies between \/2 and 2. 
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The total mass of the spectral measure on the reference simplex has a simple geometric 
interpretation. Assume that the reference norm is ii, i.e. ||x|| = Xi + ■ ■ ■ + Xd ioi x & K. If 
rj e S+, then E?7i + ■ ■ ■ + E?]^ = 1, so that the £i-norm of E 77 is 1. Since cEr] = (1, . . . , 1) 
in Theorem 12.11 we have c = d, i.e. the spectral measure has the total mass d. 

The weak convergence of simple max-stable random vectors can be interpreted as con- 
vergence of the corresponding max-zonoids. 

Theorem 4.3. Let ^, ^1, ^2? • • • be a sequence of simple max-stable random vectors with spec- 
tral measures a, ai, o"2, . . . and dependency sets K, Ki, K2, . . .. Then the following statements 
are equivalent. 

(i) ^„ converges in distribution to ^; 

(ii) o"„ converges weakly to a; 

(iii) Kn converges in the Hausdorff metric to K. 

Proof. The equivalence of (i) and (ii) is well known, see [HI Cor. 6.1.15]. 

The weak convergence of (T„, the continuity of h{Aa,x) for a G §+ and (12.41) imply that 
the support function of Kn converges pointwisely to the support function of K. Because 
dependency sets are contained inside the unit cube and so are uniformly bounded, their con- 
vergence in the Hausdorff metric is equivalent to the pointwise convergence of their support 
functions. 

The Hausdorff convergence of Kn to K implies the pointwise convergence of their support 
functions and so the pointwise convergence of the cumulative distribution functions given by 
(12. 6p . The latter entails that ^„ converges in distribution to ^. □ 

A random vector C G E belongs to the domain of attraction of a simple max-stable 
distribution if and only if the measure 

a,(A) = .P{^GA,||C||>s}, AC§+, (4.1) 

converges weakly as s — > cxo to a finite measure on §+, which then becomes the spectral 
measure of the limiting random vector, see [21 (8.95)]. The equivalence of (ii) and (iii) in 
Theorem 14.31 implies the following result. 

Proposition 4.4. A random vector ^ G E belongs to the domain of attraction of a simple 
max-stable random vector ^ with spectral measure a if and only if the max-zonoids generated 
by as from A4.1\) converge in the Hausdorff metric as s ^ 00 to the max-zonoid generated 
by cr. 
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5 Copulas and association 



The dependency structure of a distribution with fixed marginals can be explored using the 
copula function C defined on I'^ = [0, 1]'' by the following equation 

F{x) = F{x^, . . . , X,) = . . . , Frf(xd)) , 

where Fi, . . . ,Fd are the marginals of F, see [26]. In case of a simple max-stable distribution, 
we obtain 

C(mi, ...,Ud)= exp{-||(-logMi, . . . ,-\ogUd)\\K} ■ (5.1) 

Theorem 5.1. The function ( 15. Jj) is a copula function if and only if K is a niax-zonoid. 

Proof. The sufficiency is trivial, since the right-hand side of (15.11) can be used to construct a 
max-stable distribution. In the other direction, we can substitute into C the $i-distribution 
functions, i.e. Ui = e~^^^\ This yields a multivariate cumulative distribution function given 
by F{x) = e""^*"-^. The result then follows from Theorem 13. 1[ □ 

Note that ( 15. ip in the bivariate case appears in [H] . A rich family of copulas consists of 
the Archimedean copulas that in the bivariate case satisfy ip{C{xi,X2)) = <^{xi) + <^{x2) for 
a strictly decreasing continuous function and all xi,X2 G [0, 1], see [251 Ch. 4]. Using (15. ip . 
it is easy to see that in this case a;2)||_ft:) = "ipixi) + ip{x2) for a monotone increasing 

continuous function ifj and all Xi,X2 > 0. It is known (see [26| Th. 4.5.2] and [15]) that 
all Archimedean copulas that correspond to max-stable distributions are so-called Gumbel 
copulas, where ip{t) = t^. Thus, (15. ip is an Archimedean copula if and only if K is £p-ball 
with p G [1, oo], see Example | 



The bivariate copulas are closely related to several association concepts between random 
variables, see [2lj. The Spearman correlation coefficient is expressed as ps = 12 J — 3, where 

J= [ [ C{ui,U2)duidu2= [ [ e-ll(-^°s"^'-^°s"2)l'^dMidu2 



^0 ^0 ^0 

oo poo 



i / 0, oo) VII^'I'Mx , 



and 



It is possible to calculate J by changing variables x = r{t, 1 — t) with r > and t G [0, 1], 
which leads to the following known expression 

1 1 

J = - I rrr^dt = I — — —dt 



4 7o \\it,i-t)\\i J, ii + \\{t,i-t)\\Ky 
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see [19]. The following proposition is useful to provide another geometric interpretation of 
Ps and also an alternative way to compute J. 



Proposition 5.2. li L is a convex set in MJ^, then 

[ e-^(^'")cia; = r(ci+l)yd(L°), 

i[0,oo)<* 

where Vd{-) is the d-dimensional Lebesgue measure, L° is the polar set to L and T is the 
Gamma function. 

Proof. The proof follows the argument mentioned in [571 P- 2173]. Let ( be the exponentially 
distributed random variable of mean 1. Then 



C>h{L,x) 

[0,oo)<* J[0,oo)<* 

= EVd{{x G E : h{L,x) < (}) 
= EVd{Cn = Vd{L")EC''. 

It remains to note that EC"^ = r(rf+ 1). □ 
Thus, in the planar case 

PS = 3{2V2iL'^) - 1) . 

As a multivariate extension, an affine function ps = c{Vd{L°) — a) of the d-dimensional 
volume of L° may be used to define the Spearman correlation coefficient for a d-dimensional 
max-stable random vector with unit Frechet marginals. By considering the independent 
case L = K = with ps = and using the fact that the volume of the unit cross- 
polytope L° is {d\)~^, we see that ps = c{d\Vd{L°) — 1) for some constant c > 0. The 
choice c = (d + l)/(2'^ — d — 1) ensures that = 1 in the totally dependent case, where 
=27(d + l)!. 

The Kendall correlation coefficient of a bivariate copula C is given by 

•1 rl 



T = A [ [ C{Zi,Z2)dC{Zi,Z2) 

Jo Jo 



/ ——C{zi,Z2)T:—C{zi,Z2)dZidZ2. 

Jo dZ2 



see [Ml (2.3)]. By (15.11) . the partial derivatives of C can be expressed using partial derivatives 
of the support function of K. The directional derivative of the support function h{K, x) at 
point X in direction u is given by h{F{K,x),u), where 

F{K,x) = {yeK: {y,x) = h{K,x)} 
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is the support set of K in direction x, see [33l Th. 1.7.2]. Thus the partial derivatives of 
h{K,x) are given by 

^^^1^ = h{F{K, x), (1, 0)) = y.{K, x) , z = l,2, 

where yi{K, x) and y2{K, x) are respectively the maximum first and second coordinates of the 
points from F{K, x). If the dependency set K is strictly convex in (0, oo)^, i.e. the boundary 
of K inside (0, cxd)^ does not contain any segment, then F{K,x) = {{yi{K,x),y2{K,x))} is 
a singleton for all a; G E. In this case denote 

y{K, x) = yi{K, x)y2{K, x) . 

By using (15. ip and changing variables we arrive at 

r = 1 - 4 [ e~2"^ll^y(fs:,x)dx. 

J[0,oo)2 

The fact that y{K,tx) = y{K,x) and a similar argument to Proposition 15.21 yield that 



r = 1-2 / y{K,x)dx. (5.2) 
Jk° 

For instance, r = 1/2 if ^ has the logistic distribution with parameter 1/2, i.e. i| ■ is the 
Euclidean norm. By changing variables x = (t, 1 — t)r, we arrive at 



^_ [' yi{K, {t,l-t))y2{K, (t,l-t)) ^ 



which also corresponds to fT9^, Th. 3.1]. 



The Pearson correlation coefficient for the components of a bivariate simple max-stable 
random vector is not defined, since the unit Frechet marginals are not integrable. However 
it is possible to compute it for the inverted coordinates of ^. 

Proposition 5.3. If ^ is a simple max-stable bivariate random vector, then E(^f ^^2^^) = 
2V2{K°), and the covariance between and ^2"^ is 2V2{K°) - 1. 



Proof. Integrating by parts, it is easy to see that 



mi^i2^) = / / F{x*)dx,dX2 = / 



The result follows from Proposition 15.21 and the fact that E(^]^ ^) = E(,^2 ^) = 1- ^ 
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Proposition 15.31 corresponds to the formula 



for the covariance obtained in [35] for the exponential marginals. 

Extending this concept for the higher-dimensional case, we see that the covariance matrix 
of ^* is determined by the areas of polar sets to the 2-dimensional projections of K and 

dW.jK'^) - 1 
d\-l 

can be used to characterise the multivariate dependency of a simple d-dimensional max-stable 
random vector so that p varies between zero (complete independence) and 1 (complete 
dependence) . 

Example 5.4. Assume that \\x\\k = \\x\\p is the £p-norm with p > 1, i.e. the corresponding 
has the logistic distribution with parameter a = 1/p, see Example 13.31 The volume of the 
ip-hall {x eW^ : \\x\\p < 1} equals 

(2r(l + l/p)Y 



T{l + d/p) 



see [271 P- H]- Thus, the volume of K° is 2 '^Vdip) and the multivariate dependency of can 
be described by 

^_ 1 f^,{T{l + llp)Y 



d\-l\' r(l + d/p) 
If (i = 2, then p = aB{a, a) — 1 with B being the Beta-function. 

The tail dependency index for ^ = (,^i, ^2) with identical marginal distributions supported 
by the whole positive half-line is defined as 

X = lim P{6 > t|6 > t} . 

t^oo 

An easy argument shows that x = 2 — ||(1, l)||_ft' if ^ has a simple max-stable distribution 
with dependency set K, cf [4J. 

It is easy to see that ^ has all independent coordinates if and only if ||(1, . . . , 1)||_r- = d 
and the completely dependent coordinates if and only if ||(1, . . . , 1)\\k = 1, cf [31] and [21 
p. 266]. It is well known [2, p. 266] that the pairwise independence of the coordinates of ^ 
implies the joint independence. Indeed, the spectral measure of the set m G S_|_ such that at 
least two coordinates of u are positive is less than the sum of a{u G S+ : Ui > 0, Uj > 0} 
over all i ^ j- Each of these summands vanishes, since 



Xi + Xj — / max {ukXk)(y{du) = / {{uiXi + ujXj) — {uiXi V UjXj))a{du) = 

by the pairwise independence, where x has all vanishing coordinates apart from Xi and Xj. 
This leads to the following property of max-zonoids. 
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Proposition 5.5. If K is a inax-zonoid with all its two-dimensional projections being unit 
squares, then K is necessarily the unit cube. 

6 Complete alternation and extremal coefficients 

Consider a numerical function / defined on a semigroup S with a commutative binary op- 
eration +. For n > 1 and Xi, . . . ,Xn G S define the following successive differences 

^xJix) = f{x) - f{x + Xi) , 

• • ■ AxJ{x) = A^„_j • • ■ A^J{x) - A^„_, ■ ■ ■ A^J{x + Xn) ■ 

The function / is said to be completely alternating (resp. monotone) if all these successive 
difference are non-positive (resp. non- negative) , see [U Sec. 4.6] and [23^ Sec. 1.1.2]. We will 
use these definitions in the following cases: S is the family of closed subsets of M'^ with the 
union operation, S is M.'^ or E with coordinatewise minimum or coordinatewise maximum 
operation. Then we say shortly that the function is max-completely alternating or monotone 
(resp. min-completely or union-completely). 

Every cumulative distribution function F is min-completely monotone. This is easily 
seen by considering the random set X = {^} where ^ has the distribution F, then noticing 
that F{x) = P{X n Lj. = 0} = Q^L^) with being the complement to x + {—oo, OY is the 
avoiding functional of X, and finally using the fact that 

F(min(x, y)) = P{X n (L, U Ly) = 0} = Q(L, U Ly) 

together with the union-complete monotonicity of Q, see [23] Sec. 1.6]. 

Theorem 6.1. A convex set K G K is a max-zonoid if and only if h{K,x) is a max- 
completely alternating function of x . 

Proof. It follows from [21 Prop. 4.6.10] that a function f{x) on a general semigroup is com- 
pletely alternating if and only if F{x) = e"*^^"^'-* is completely monotone for all t > 0. Since 
{mm{x,y)y = max{x*,y*) for x, |/ G E, the function x i— > f{x*) is max-completely alternat- 
ing on E if and only if x i— *• f{x) is min-completely alternating. 

If K is the max-zonoid corresponding to a simple max-stable random vector with distri- 
bution function F{x) = e~^^^'^*\ then F* is also a cumulative distribution function (and so 
is min-completely monotone) for each t > 0. Thus, h{K,x*) is min-completely alternating, 
whence h{K, x) is max-completely alternating. 

In the other direction, if h is max-completely alternating, then F(x) = q-K^''-^*) is min- 
completely monotone, whence it is a cumulative distribution function. The corresponding 
law is necessarily semi-simple max-stable, so that K is indeed a max-zonoid. □ 
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Theorem 16.11 can be compared with a similar characterisation of classical zonoids, where 
the complete alternation of h[K, x) and monotonicity of e~'^^^'^^ are understood with respect 
to the vector addition on M.'^, see p. 194]. 

The remainder of this section concerns extensions for the support function defined on a 
finite subset of E. 

Theorem 6.2. Let M be a finite set in E, which is closed with respect to coordinatewise 
maxima, i.e. uM v & M for all u,v ^ M. Assume that for each u,v & M, we have tu < v if 
and only if u <v and t <1. Then a non-negative function h on M can be extended to the 
support function of a max-zonoid if and only if h is max- completely alternating on M. 

Proof. The necessity trivially follows from Theorem 16. 1[ To prove the sufficiency we ex- 
plicitly construct (following the ideas of [3T]) a max-stable random vector ^ such that the 
corresponding norm coincides with the values of h on the points from M. 

For any set A C M, let \/A denote the coordinatewise maximum of A. Furthermore, 
define T{A) = h{\/ A) / hiy M) . Since h is max-completely alternating, T is union-completely 
alternating on subsets of M. The Choquet theorem [231 Th. 1.1.13] implies that a union- 
completely alternating function on a discrete set is the capacity functional T{A) = P{Xr\A ^ 
0} of a random closed set X C M. Define c„ = /i(VM)P{VX = u} for uE M. 

Let (u, u E M, be the family of i.i.d. unit Frechet random variables which are also chosen 
to be independent of X and let ^ be the coordinatewise maximum of Cuu(u over all u G M. 
It remains to show that ^ has the required distribution. Consider an arbitrary point v G M. 
By the condition on M, tu < v is possible for some t > if and only ii u < v and t < 1. 
Thus, 

p{e < ^} = n < = n < 1} = exp i - 

ueM u£M, u<v \ uGM, u<v 

= exp {-/i(VM)P{X n{u: u < v} ^ ^}} 
= exp {-h{yM)T{{u : u < v})} = exp {~h{v)} . 

□ 

A simple example of set M from Theorem 16.21 is the smallest set which contains all basis 
vectors in R'^ and is closed with respect to coordinatewise maxima. Then M consists of 
the vertices of the unit cube I'^ without the origin and the values of h on M become the 
extremal coefficients. The extremal coefficients 9 a of a simple max-stable random vector 
^ = (^1, . . . , ^d) are defined from the equations 

P{maxe, <z} = (P{ei < z])'^ , ^ > , A C {1, . . . , 4 , (6.1) 

see [31], [32]. Since the marginals are unit Frechet, it suffices to use (16. ip for z = 1 only. If 
= YjieA'^i^ then ([3?I]) implies 

9a = KK.ca) = \\eA\\K- 
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Every nonempty set A C {1,. . . ,d} can be associated with the unique vertex of the unit 
cube I'^ \ {0}. The consistency condition for the extremal coefficients follows directly from 
Theorem 16.21 and can be formulated as follows. 

Corollary 6.3. A family of non-negative numbers 9 a, Ad {1, . . . , d}, is a set of extremal 
coefficients for a simple max-stable distribution if and only if 6'0 = and 9a is a union- 
completely alternating function of A. 

This consistency result for the extremal coefficients has been formulated in [31] as a set 
of inequalities that, in fact, mean the complete alternation property of 9a- 

7 Operations with dependency sets 

Rescaling. For a dependency set K and Ai, . . . , > define 

\K = {(AiXi, . . . , XdXd) : X = (xi, ...,Xd)eK}. (7.1) 
Then g-'^C^^.^*) is the cumulative distribution function of A*^ = (^i/Ai, . . . ,^rf/Ad). 

Projection. If ^' denotes the vector composed from the first /^-coordinates of d-dimensional 
vector ^ with the dependency set K, then 

P{^' < (xi, Xk)} = exp{-||(xi, . . . , Xfc, cx), . . . , oo)*||k} 
= exp{-||(xi, . . .,Xky\\K'} , 

where K' is the projection of K onto the subspace spanned by the first k coordinates in 
M'^. Thus, taking a sub-vector of C, corresponds to projecting of K onto the corresponding 
coordinate subset. Recall Proposition 15.51 which says that if all two-dimensional projections 
of K are squares, then K is necessarily the cube. 

Proposition 7.1. If h is the subspace spanned by some coordinate axes in M.'^, the projection 
of K onto L coincides with fl L. 

Proof. By definition, K = cEZ\^. Then it suffices to note that the projection of zl^ on L 
equals fl L. Indeed every selection of fl L can be associated with the projection of a 
selection of Z^^. □ 

An interesting open question concerns a reconstruction of K from its lower-dimensional 
projections. In various forms this question was discussed in [221 Sec. 3.5.6] and [21], Sec. 4.7]. 

Cartesian product. If K' and K" are two dependency sets of simple max-stable random 
vectors ^' and C," with dimensions d' and d" respectively, then the Cartesian product K' x 
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K" is the dependency set corresponding to the max-stable random vector ^ obtained by 
concatenating of independent copies of ^' and ^" . Indeed, if x = (x', x"), then 



Minkowski sum. If K' and K" are dependency sets of two independent max-stable ran- 
dom vectors ^' and ^" of dimension d, then the weighted Minkowski sum K = XK'+{1 — X)K" 
with A G [0, 1] is the dependency set of the max-stable random vector 



It is possible to generalise the Minkowski summation scheme for multivariate weights. 
Consider K = \K' + (1 — X)K" for some A G [0, l]*^, where the products of vectors and sets 
are defined in (17.11) . Then \\x\\k = ||A*x||x/ + ||(1 — \)*x\\k", so that (17.21) also holds with 
the products defined coordinatewisely. 

Example 7.2. If and ,^2 are independent with unit Frechet distributions and ai, 0^2 G [0, 1], 
then setting A = (ai, 1 — 02) we obtain the max-stable random vector 

i = {aiii V (1 - ai)6, (1 - "2)^1 V ^26) 

with the dependency set K = conv{(0, 0), (0, 1), (1, 0), (ai, 1), (1, 02)}. If ai = 02 = «, then 
^ has the Marshall-Olkin distribution, cf [121 Ex. 4.1.1]. 

Example 7.3 (Matrix weights). Let aij, i = 1, . . . ,m, j = 1, . . . ,d, be a matrix of positive 
numbers such that J^iLi^^ij = 1 all j. Furthermore, let d, . . . Xm be i.i.d. random 
variables with $i-distribution. Define ^ = (^1, . . . ,^d) by 



P{^ <x} = exp{-h{K' X K",x)} = exp{-h{K',x') - h{K'\x")} 
= P{C < x'}P{C < x"}. 



e = (Anv((i-A)n- 



(7.2) 



The cumulative distribution functions of ^' 



and ^ are related as 



F^{x) = F^,{x)''F^.{xY^'^K 



max (ittij , j = 1, . . . 



cf fl2[ Lemma 4.1.2]. Then ^ is simple max-stable with the corresponding norm 



m 




i.e. its dependency set is = Z\(, 



) + ■■■ + A< 
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Power sums. A power-mean of two convex compact sets K' and K" containing the origin 
in their interior is defined to be a convex set K such that 

h{K, xY = Xh{K', xY + (1 - X)h{K, xf , (7.3) 

where A G [0, 1] and p > 1, see |T1]. The power-mean definition is apphcable also if K' and 
K" satisfy fl2.8p . despite the fact that the origin is not their interior point. In the plane, 
the power sum is a dependency set if K' and K" satisfy fl2.8p . Therefore, the power sum of 
dependency sets leads to a new operation with distributions of bivariate max-stable random 
vectors. For instance, if K' is the unit cross-polytope and K" is the unit square, then, for 
p = 2, 

\\x\\k = {{xi + X2 f + (max(xi, X2))^)^''^ . 

Minkowski difference. Let K' and K" be two dependency sets. For any A > define 

L = K' - \K" = {x: x + XK" C K'} . 

If the spectral measures a' and a" of K' and K" are such that a = a' — Xa" is a non-negative 
measure, then L is a max-zonoid regardless of the dimension of the space. The negative 
logistic distribution from Example 13.51 illustrates this construction. 

Convex hull and intersection. In the space of a dimension d > 3 the convex hull or 
intersection of dependency sets do not necessarily remain dependency sets. However, on the 
plane this is always the case. 

Let K' and K" be the dependency sets of bivariate simple max-stable random vectors 
C and f . Since h{coiav{K' U K"),x) = h{K',x) V h{K",x), the dependency set K = 
conv(ir' U K") corresponds to a max-stable random vector ^ such that 

P{^ <x} = min(P{^' < x}, P{C < x}) , xe [0, oo)^ . 

The intersection of two planar dependency sets also remains the dependency set and 
so yields another new operation with distributions of simple max-stable bivariate random 
vectors. 

Duality. If the polar to the dependency set ii' of ^ is a max-zonoid, then the corresponding 
simple max-stable random vector ^° is said to be the dual to ^. In the plane, the polar to a 
max-zonoid is max-zonoid; it is not known when it holds in higher dimensions. This duality 
operation is a new operation with distributions of bivariate max-stable random vectors, see 
also Example 13.21 

8 Infinite dimensional case 

It is possible to define the dependency set for max-stable stochastic processes studied in 
[3 [121 [IS!- The spectral representation Prop. 3.2] of a sample continuous max-stable 
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process ^(t), t E S, on a compact metric space S with unit Frechet marginals yields that 

-logP{e</}= / \\9/f\\ooda{g), 

where S+ is the family of non-negative continuous functions g on S that their maximum value 
\\g\\oo equals 1, and a is a finite Borel measure on S+ such that Jg^ gda{g) is the function 
identically equal to 1. 

The corresponding dependency set is the set in the space of finite measures with the 
total variation distance, which is the dual space to the family of non-negative continuous 
functions. For a continuous function g, define Ag to be the closed convex hull of the family 
of atomic measures g{x)5x for x G 5. Then the dependency set is the expectation of cA^, 
where c = cr(S+) and rj is distributed according to the normalised a. 
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